Lepton flavor violating radion decays in the 
Randall- Sundrum scenario 



E. O. Iltan, *B. Korutlu f 

Physics Department, Middle East Technical University 

Ankara, Turkey 



Abstract 

We predict the branching ratios of the lepton flavor violating radion decays r — > e fj, , 
r — ► e^r 1 * 1 and r — > in the two Higgs doublet model, in the framework of the 

Randall- Sundrum scenario. We observe that their branching ratios are at most of the 
order of 10 -8 , for the small values of radion mass and they decrease with the increasing 
values of m r . Among these processes, the r — > decay would be the most suitable 

one to measure its branching ratio. 
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1 Introduction 



The hierarchy problem between weak and Planck scales could be explained by introducing the 
extra dimensions. One of the possibility is to pull down the Planck scale to TeV range by 
considering the compactified extra dimensions of large size [T] . The assumption that the extra 
dimensions are at the order of submilimeter distance, for two extra dimensions, the hierarchy 
problem in the fundamental scales could be solved and the true scale of quantum gravity would 
be no more the Planck scale but it is of the order of electroweak (EW) scale. This is the 
case that the gravity spreads over all the volume including the extra dimensions, however, the 
matter fields are restricted in four dimensions, so called four dimensional (4D) brane. Another 
possibility, which is based on the non-factorizable geometry, is introduced by Randall and 
Sundrum [2] (the RSI model) and, in this scenario, the extra dimension is compactified to 
S 1 /Z2 orbifold with two 4D brane boundaries. Here, the gravity is localized in one of the 
boundary, so called the Planck (hidden) brane, which is away from another boundary, the TeV 
(visible) brane where we live. The size of extra dimension is related to the vacuum expectation 
of a scalar field and its fluctuation over the expectation value is called the radion field (see 
section 2 for details). The radion in the RSI model has been studied in several works in the 
literature [3]-[l3] (see [Ej for extensive discussion). 

In the present work, we study the possible lepton flavor violating (LFV) decays of the radion 
field r. The LFV interactions exist at least in one loop level in the extended standard model 
(SM), so called ^SM, which is constructed by taking neutrinos massive and by permitting the 
lepton mixing mechanism [15J. Their negligibly small branching ratios (BRs) stimulate one to 
go beyond and they are worthwhile to examine since they open a window to test new models 
and to ensure considerable information about the restrictions of the free parameters, with the 
help of the possible accurate measurements. The LFV interactions are carried by the flavor 
changing neutral currents (FCNCs) and in the SM with extended Higgs sector (the multi Higgs 
doublet model) they can exist at tree level. Among multi Higgs doublet models, the two Higgs 
doublet model (2HDM) is a candidate for the lepton flavor violation. In this model, the lepton 
flavor violation is driven by the new scalar Higgs bosons S, scalar h° and pseudo scalar A , and 
it is controlled by the Yukawa couplings appearing in lepton-lepton-S vertices. 

Here, we predict the BRs of the LFV r decays in the 2HDM, in the framework of the RSI 
scenario Q- The BRs of these processes are sensitive to the new Yukawa couplings arising with 

1 The flavor conserving decays of the radion can exist in the tree level and the theoretical values of the BRs 
are free from the loop suppressions. For the heavy flavor output rr the BR reaches to the numerical values of 
the order of 10~ 2 for the small values of the radion mass, m r < 200 GeV. In the case of light flavor output 
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the addition of the new Higgs doublet and they enhance with the increasing values of the these 
couplings. Therefore, besides their important role in the construction of the FCNC in the tree 
level and the existence of the LFV interactions, these Yukawa couplings play a crucial role in 
the enhancement of the theoretical values of the BRs. In our calculations, we observe that 
the BRs of the processes we study are at most of the order of 10~ 8 , for the small values of 
radion mass m r and for the large values of the Yukawa couplings, and, their sensitivities to m r 
decrease with the increasing values of m r . Among the LFV decays we study, the r — > ^ 
decay would be the most suitable one to measure its BR. 

The paper is organized as follows: In Section 2, we present the effective vertex and the 
BRs of LFV r decays in the 2HDM, by respecting the RSI scenario. Section 3 is devoted 
to a discussion and to our conclusions. In the appendix, we present the interaction vertices 
appearing in the calculations. 

2 The LFV RSI radion decay in the 2HDM 

The RSI model is an interesting candidate in order to explain the well known hierarchy problem. 
It is formulated as two 4D surfaces (branes) in 5D world in which the extra dimension is 
compactified into S 1 /Z2 orbifold. In this model, the SM fields are assumed to live on one of the 
brane, so called the TeV brane. On the other hand, the gravity peaks near the other brane, so 
called the Planck brane and extends into the bulk with varying strength. Here, 5D cosmological 
constant is non vanishing and both branes have equal and opposite tensions so that the low 
energy effective theory has flat 4D spacetime. The metric of such 5D world reads 

ds 2 = e- 2A{y) ri^dx»dx u -dy 2 , (1) 

where A(y) = k \y\, k is the bulk curvature constant, y is the extra dimension parametrized as 
y = R9. The exponential factor e~ kL with L = Rn, is the warp factor which causes that all 
the mass terms are rescaled in the TeV brane. With a rough estimate L ~ 30/ k, all mass terms 
are brought down to the TeV scale. The size L of extra dimension is related to the vacuum 
expectation of the field L(x) and its fluctuation over the expectation value is called the radion 
field r. In order to avoid the violation of equivalence principle, L(x) should acquire a mass 
and, to stabilize r, a mechanism was proposed by Goldberger and Wise [3], by introducing a 
potential for L(x). Finally, the metric in 5D is defined as [I] 

ds 2 = e -2My)-2F(x) Vfiudx » dx » _(i + 2 F(x))dy 2 , (2) 

yy (ee) the numerical value of the BR is of the order of 10~ 4 (10 -9 ) for m r < 200 GeV (see [9] for the explicit 
expressions of these decay widths) . 
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where the radial fluctuations are carried by the scalar field F(x) 

1 



F(x) 



-kL 



r(x) 



(3) 



V6M Pl e 

Here the field r(x) is the normalized radion field (see [5]). At the orbifold point 9 = n (TeV 
brane) the induced metric reads, 



ind 



e - 2 A(L)- 2 J,(,)^ 



(4) 



Here the parameter 7 reads 7 



with A = Mpie kL and v is the vacuum expectation 



value of the SM Higgs boson. The radion is the additional degree of freedom of the 4D effective 
theory and we study the possible LFV decays of this field. 

The FCNCs at tree level can exist in the 2HDM and they induce the flavor violating (FV) 
interactions with large BRs. The FV r decays, r — > can exist at least in one loop level in 
the framework of the 2HDM. The part of action which carries the interaction, responsible for 
the LFV processes reads 



S 



Y 



d^xyj -g^ir^kL^EjR + ^U L ^E jR + h.c. 



(5) 



where L and R denote chiral projections L(R) = 1/2(1 =F 75), 0j for i = 1,2, are two scalar 
doublets, In (Ejn) are lepton doublets (singlets), ifjU and rj^, with family indices i,j , are 
the Yukawa couplings and ^ induce the FV interactions in the leptonic sector. Here g 



ind 

is the determinant of the induced metric on the TeV brane where the 2HDM particles live. 
Here, we assume that the Higgs doublet <pi has a non-zero vacuum expectation value to ensure 
the ordinary masses of the gauge fields and the fermions, however, the second doublet has 
no vacuum expectation value, namely, we choose the doublets <pi and 02 and their vacuum 
expectation values as 



1 

V2 







and 



< 



( V2x 



wo 
V2 V v 



y/2 \ Hi + iH 2 



(6) 



< 02 >= . 



(7) 



This choice ensures that the mixing between neutral scalar Higgs bosons is switched off and it 
would be possible to separate the particle spectrum so that the SM particles are collected in the 



In the following, we replace £ B with £jy where "N" denotes the word "neutral". 
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first doublet and the new particles in the second one[]. The action in eq. flS]) is responsible for 
the tree level S — h — h (h and I2 are different flavors of charged leptons, S denotes the neutral 
new Higgs boson, S = h°, A ) interaction (see Fig. 1-d, e) and the four point r — S — l\ — 1% 
interaction (see Fig. 1-c) where r is the radion field. The latter interaction is coming from 



the determinant factor y—g md = e~ 4j4(X - ) ~ 4 ^ t<kX \ Notice that the term e~ 4j4 ^^ in y—g' md is 
embedded into the redefinitions of the fields on the TeV brane, namely, they are warped as 
S — > e A ^ S warp , I — ► e 3A( * L ^ 2 l W arp and in the following we use warped fields without the warp 
index. 

On the other hand, the part of new scalar action 

S 2 = J d 4 xJ^T d (g ind ^(D, (j>^D v( j> 2 -m 2 s 4fo^J (8) 

leads to 



S' 2 = I J d A x je" 2 ? r V ^ (d, h°d u h° + d, A°d u A ) - e" 4 ? 



m 2 h° h° + m 2 A0 A A ) \ , (9) 



which carries the S — S — r interactional (see Fig. 1-b). 

Finally, the interaction of leptons with the radion field is carried by the action (see [6]) 

S3 = Jd i x^^ d ^g md ^h^DJ-m l ll ] j, (11) 

where 

D^ = d,i + ^ w ; b j: ab i, (12) 

with S afe = \ [7 a , jb\. Here w ab is the spin connection and, by using the vierbein fields e* it can 
be calculated (linear in r) as 

w f = -ld„r(e ub e*-e va el) . (13) 



3 Here Hi (H2) is the well known mass eigenstate h° (A ). 

4 In general, there is no symmetry which forbids the curvature-scalar interaction, 

5 e = J d?x^-g ind tK Hi H , (10) 

where £ is a restricted positive parameter and H is the Higgs scalar field [T3l [14] . This interaction results in 
the radion- (SM or new) Higgs mixing which can bring a sizeable contribution to the physical quantities studied. 
Here, we assume that there is no mixing between first and second doublet and only the first Higgs doublet has 
vacuum expectation value. Therefore, we choose that there exists a mixing between the radion and the SM 
Higgs field, but not between the radion and the new Higgs fields. This is the case that the lepton flavor violation 
is not affected by the mixing since the SM Higgs field is not responsible for the FCNC current at tree level. 
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Notice that the vierbein fields are the square root of the metric and they satisfy the relation 

e» e av = g ind ^ . (14) 

Using eqs. (fTIl) - (IT4"j) . one gets the part of the action which describes the tree level I — I — r 
interaction (see Fig.l-a) as 

S' 3 = (d A x\-?>-rliBl-?>—li8rl + A-mirTl\. (15) 
J { v 2v v J 

At this stage, we would like to discuss briefly the case that the gauge fields and the leptons 

are also accessible to the extra dimension. With the addition of Dirac mass term moill) = 

m^-j^t (A'(y) = dA ^ ) to the lagrangian of bulk fermions, the fermion mass hierarchy can 

be obtained since they have different locations, regulated by a localization parameter c = 

m/k, in the extra dimension (see for example [THl d3 [TSl ) . Furthermore, the SM fermions 

are obtained by considering SU{2) L doublet ipL and singlet ipR with separate Z 2 projection 

conditions: Z 2 i/jr = I^r and Z 2 %1)l = —Ih^L (see for example [20J) and the localization 

parameters cl and cr assigned to each left and right handed fermion. If cl-r > | {cl-r < f), 

the zero mode left;right handed fermions are localized near the hidden (visible) brane and, for 

C L,R = h they have constant profiles. Notice that the parameters cl and cr should be restricted 

such that the observed masses of fermions are obtained. In addition to the SM fermions (the 

zero mode ones) the KK modes of left and right handed fermions arise and these modes depend 

on the localization parameters and cr. In this scenario, the BRs of the LFV radion decays 

become sensitive to the different locations of left and right handed (zero mode and KK mode) 

leptons. This sensitivity is due to the couplings coming from r — ^(i?) — ^t?(l)' r — ^x(R) — ^x(r) 

and S — ^2(_r) ~ ^iR(L) interactions, where S = h°,A and n = 0,1,2,.... Here we can take 

two different possibilities: The first possibility is that the existence of the LFV is based on 

the different locations of different flavors (see for example [2~Tl I2"2"l [23] )• In this case the new 

Yukawa couplings, appearing with the interaction of zero mode leptons and new Higgs bosons, 

sensitive to the locations of the corresponding leptons and the sources of the flavor violation are 

explicitly the localization parameters. If the leptons are far from (near to) the visible brane the 

couplings are suppressed (enhanced). Furthermore, the strengths of the couplings of KK mode 

leptons and new Higgs bosons are regulated by the locations of the lepton fields. In the second 

possibility, the flavor violation is carried by the new Yukawa couplings in four dimensions by 

assuming that the localization effects are embedded into the definition of new Yukawa couplings 

and, with the additional effects coming from the KK mode leptons, the role of extra dimension 

becomes the enhancement in the physical quantities . In both scenarios, the couplings of the 
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radion with the leptons are regulated by the localization of the zero mode leptons and there 
exists an additional effect coming from the radion-KK mode lepton vertices. At first sight, it 
is expected that if the zero leptons are far from (near to) the visible brane these couplings are 
suppressed (enhanced). However, one needs a detailed analysis to take into account the bulk 
lepton contributions and their localization effects. 

Now, we are ready to calculate the matrix element for the LFV radion decay. The decay 
of the radion to leptons with different flavors exits at least in one loop order, with the help 
of internal new Higgs bosons S = h°, A . The possible vertex and self energy diagrams are 
presented in Fig. [21 After addition of all these diagrams, the divergences which occur in the 
loop integrals are eliminated and the matrix element square for this decay is obtained as 

\M\ 2 = 2(m 2 r -(m lI+ m lt ) 2 )\A\\ (16) 

where 

A fself | rself | rvert i rvert i rvert i rvert f1'7\ 

^ — ho + J AO + Jh° + J A + Jh°h + J A A ) [i<) 



and the explicit expressions of the functions appearing in eq. (fTTj) are given in Appendix B. 
Notice that in eq. (T231) the flavor changing couplings Cn^i represent the effective interaction 
between the internal lepton i, {i = e, /x, r) and the outgoing j = 1 (j = 2) lepton (anti lepton). 
Here, we choose the couplings £§j i real. 

Finally, the BR for r — > l± 1% can be obtained by using the matrix element square as 

1 \M\ 2 



BR(r -> /J" 1' 2 



lQ-Knir T r 



where T r is the total decay width of radion r. In our numerical analysis, we consider the BR 
due to the production of sum of charged states, namely 

+ +x T(r -> (kh + l 2 h)) , . 

BR(r^lflf) = ^ V ;; . (19) 

i r 

3 Discussion 

In four dimensions, the higher dimensional gravity is observed as it has new states with spin 
2,1 and 0, so called, the graviton, the gravivector, the graviscalar. These states interact with 
the particles in the underlying theory. In the RSI model with one extra dimension, the spin 
gravity particle radion r interacts with the particles of the theory (2HDM in our case) on the 
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TeV brane and this interaction occurs over the trace of the energy-momentum tensor with 
the strength 1/A r , 



where A r is at the order of TeV. The radion interacts with gluon (g) pair or photon (7) pair in 
one loop order from the trace anomaly. For the radion mass m r < 150 GeV, the decay width 
is dominated by r — > gg. For the masses which are beyond the WW and ZZ thresholds, the 
main decay mode is r — ► WW. In the present work, we study the possible LFV decays of the 
RSI radion in the 2HDM and estimate the BRs of these decays for different values of radion 
masses. We take the total decay width T r of the radion by considering the dominant decays 
r — ► gg (77, //, W + W~ , ZZ, SS) where S are the neutral Higgs particles (see [9] for the explicit 
expressions of these decay widths). Here, we include the possible processes in the r r according 
to the mass of the radion. 

The flavor violating r decays r — > can exist at least in one loop level, in the framework 
of the 2HDM and the flavor violation is carried by the Yukawa couplings £jvijl In the version 
of 2HDM where the FCNC are permitted, these couplings are free parameters which should be 
restricted by using the present and forthcoming experiments. At first, we assume that these 
couplings are symmetric with respect to the flavor indices i and j. Furthermore, we take that 
the couplings which contain at least one r index are dominant and we choose a broad range 
for these couplings, by respecting the upper limit prediction of £,§ tTfl (see [24J and references 
therein), which is obtained by using the experimental uncertainty, 10~ 9 , in the measurement of 
the muon anomalous magnetic moment and by assuming that the new physics effects can not 
exceed this uncertainty. For the coupling £§ Te , the restriction is estimated by using this upper 
limit and the experimental upper bound of BR of /1 — ► ej decay, BR < 1.2 x 10 -11 [25]. Finally, 
this coupling is taken in the range 1CT 3 — 1CT 1 GeV (see [26]). For the Yukawa coupling ^§ TT , 
we have no explicit restriction region and we use the numerical values which are greater than 
£§ TlV Throughout our calculations we use the input values given in Table ([T]). 

In Figj3]we present m r dependence of the BR (r — > r fx ). The solid-dashed lines represent 



the BR (r - r^) for ^ TT = 100 GeV, ^ = 10 GeV- £f jTT = 10 GeV, = 1 GeV. 



It is observed that the BR (r — > fj,^) is of the order of the magnitude of 10~ 8 for the large 
values of the couplings and the radion mass values ~ 200 GeV. For the heavy masses of the 
radion the BR is stabilized to the values of the order of 10 -9 . 



int 



A r M 



(20) 



5 The dimensionfull Yukawa couplings £j| -■ are denned as £ 



■E 

N,ij 
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Parameter 


Value 


"V 
m T 
m h o 

Gp 


U. 1UU y\jrc V ) 

1.78 (GeV) 
100 (GeV) 
200 (GeV) 
1.1663710" 5 (GW- 2 ) 



Table 1: The values of the input parameters used in the numerical calculations. 

FigJH is devoted to m r dependence of the BR (r — > re) and BR (r — > /i e ). The 
solid-dashed lines represent the BR (r -»■ r± e ± ) for ^ TT = 100 GeV, £f )Te = 0.1 GeV- £§ >TT = 
10 GeV, £jf re = 0.1 GeV. The small dashed line represents the BR (r — > /i ± e ± ) for ^ r „ = 
lGeV, ifi Te = 0.1 GeV. This figure shows that the BR (r — > y^) is of the order of the 
magnitude of 10~ 12 for the large values of the couplings and the radion mass values ~ 200 GeV . 
For the heavy masses of the radion, this BR reaches to the values less than 10 -14 . The BR 
(r — > /i ± e ± ) is of the order of 10~ 15 for m r ~ 200 GeV and for the intermediate values of 
Yukawa couplings. These BRs, especially BR (r — ► /i ± e ± ), are negligibly small. 

Now, we present the Yukawa coupling dependencies of the BRs of the decays under consid- 
eration, for different radion masses . 

Fig J5] represents the Cntt dependence of the BR (r — >• r^ 1 /j^) for ^ T „ = 10 GeV. The solid- 
dashed-small dashed lines represent the BR for the radion masses m r = 200 — 500 — 1000 GeV. 
This figure shows that the BR is sensitive to the radion mass and, obviously, it is enhanced two 
orders of magnitude in the range 10 GeV <Cntt— 100 GeV. 

In FigJHl we present the £§ TT dependence of the BR (r — ► e^ for £§ Te = 0.1 GeV. 
The solid-dashed-small dashed lines represent the BR for the radion masses m r = 200 — 500 — 
1000 GeV. Similar to the r — > r* ^ decay, the BR is strongly sensitive to the radion mass. 

4 Conclusion 

The LFV decays of the radion in the RSI model strongly depend on the radion mass and the 
Yukawa couplings. The BR for r — ► r fi decay is of the order of 10~ 8 for the small values of 
radion mass m r and it decreases with the increasing values of m r . On the other hand, the BRs 
for r — > r ± e (r — > / u ± e ) decays are of the order of 10~ 12 (10~ 15 ) for the small values of m r . 
These results show that, among these processes, the LFV r — > r^ ^ decay would be the most 
appropriate one to measure its BR. The most probable production of radion is due to the gluon 
fusion, gg —>■ r [9j) and, based on the existing cross section values for the radion production via 
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gluon fusion and the typical integrated luminosities expected at LHC, a rough estimation of the 
number of events for the ^ final state is missing at present. The goal being to check that 
this number of events is significant statistically, and hopefully, the forthcoming experimental 
results of this decay would be useful in order to test the possible signals coming from the extra 
dimensions and new physics which results in flavor violation. 

A The vertices appearing in the present work 

In this section we present the vertices appearing in our calculations. Here S denotes the new 
neutral Higgs bosons h° and A . 



9 



li 



Pi 



P2 



-i-f 
v 



{$1+ ^2) - 4 m k 



Pi 

P2 



^(Pi.p 2 -m|) 



(5 = A°) 



—47 

2^1) 



li 



[($ + $*) + - £jt*h 5 ; 



1 

2\/2 



(eg-ef) + (C? + ef)7 5 ' 



Figure 1: The vertices used in the present work. 
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B The explicit expressions appearing in the text 



nself 

Jh° 



fself 
J AO 



rut 
Jh° 



■vert 



vert 



fee. 
JA° 



h° 



The explicit expressions of functions f^oLoy fh° r (A°)i fh° T h° anc ^ in ec L-([I3) are given by 

1 f J^ se ^f 2 

= 128v7r^(w' h -w h ) So dX mh ° { $ (X " 1} Wh ~ V ' Zlh ^ ( 3w ' h ~ 5 Wh ) ln 1,h ° ^ 
+ (riY (x -l)w' h - nf z ih J (5 w' h - 3 w h ) In 



t self 2 
^2,11° m h° 



7 



128 v tc 2 (w' A — wa) Jo 



dx , m A o | (pY (x - 1) w A + vt z iA) (3 w A " ' " 1 ' 



5wa) In 



+ (j]Y (x -l)w' A + nf z iA ) [5w A -3 w A ) 



In 



jself 2 

M 2 



i r dx r x dy j m 

128 v vt 2 Jo Jo 1 L«f 



^ 3* (y (! - y) «4 + x 2 (4 y - 1) w h 
+ x((l-3y)w h + y(4y-3)w' h )) +5z 2 ((2y-l)w' h + (2x-l)w h ) 
3(x + y- 1) (x {Ax -3) w£ + y(4y-3) 

3 iu h «4 (x + y - 1) ((1 - y + x (4 y - 2)) w h + (1 - 2 y + x (4y - 1)) w' h 



+ 3 



(x + y - 1) ((2 x - 1) w fe + (2 y - 1) ^) 



+ Vt z ih (^x + y-l)(-A + 2w' h w h + w' 2 (8y - 3) + to; 



(8x-3)) 



- (84 + ^\((8y-3)x-3y)) 



m^o In ■ 

7 

128 W7T 2 Jo 



ver ^2 
J h° m h° 



9 r)Y (w' h {2y-l)+w h (2x-lj)-8 r)+ z ih 



dx 



l-x 



dy 



m>AO 

j ver 
. ^A° 



vY \3 4a (y (! v) w 'a ~\~ x 2 (Ay 1) w A 

+ x((l-3y)^ + y(4y-3)^)) +5^ 2 ((2y-l)^ + (2x-l) WA ) 

- 3 (x + y - 1) (x (4 x - 3) w\ + y (4 y - 3) iuJ) 

- 3w A u4(x + y-l) ((1 -y + z(4y-2))w A + (1 -2y + x(4y- l))w' A 
+ 3(x + y-l) ((2x-l)io A + (2y-l)<)j 

+ ^i" ((x + y - 1) ( - 4 + 2w' A w A + w'2 (8y - 3) + w 2 A (8 x - 3)) 

~ (8z 2 A + z 2 rA ((8y-3)x-3y)) 



n%AO In ■ 



ver m 2 
J A° A 



9 V Y(w' A (2y-l)+w A (2x-l))+8ntz iA 
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rA = cfa/ o -^{^-[^(y-l + *(l-4y))(x«, fc + y«/ h ) 

+ y(x + y-l)w' h ((Ax-l)w 2 + (Ay-l)w' 2 ) 



+ w 3 h x(x + y - 1) (Ax - 1) + (x + y - 1) (2yw' h + xw h (2 + u>; 2 (4y - 1)) )^ 

+ r ? +^(x + y-l)^((4y-l)^ 2 + (4a;-l)^ + 2) -^z r \((4y-l)a;-y + l) 

- m,o /n (i _ 6 y) + Wfc (1 - 6 x)) - 4 ^ ^ J L , 



+ y(x + y-l)<((4x-l)< + (4y-l)u/f) 



+ ^ 



x (x + y - 1) (4x - 1) + (x + y - 1) (2yw' A + xw A (2 + w 2 A (Ay- 1)) fj 
+ ^ + ^ + y - 1) ^ ((4y - 1) w'l + (Ax-l)w 2 A + 2) - z iA z 2 A ((Ay - 1) x - y + lfj 
- m A o In L ™*<° [r)?(w' A (l-6y)+w A (l-6x))+A4 z iA ^j } , (21) 



where 



^l^o (A") ~~ 1 + X 2 W^A) + x ( z ih(iA) ~ w h(A) " 1) > 

r self I | 2 '2 , / 2 '2 i\ 

^2,fc°(A°) — 1 + x ' U, ft(A) + X l Z j/i(iA) — W h(A) ~ l ) > 

l )!p(a°) = x2 w ma) + (y - 1) ( w l 2 (A) y - !) + & (y w 'h(A) + (y - 1) w ma) - y 4>m) - 1) , 

^W(A° AO) = X 2 W 2 h(A) + (1 + ^ 2 (A) (y - !)) y + ^ (! + ™l(A) (V ~ 1) + W MA) 2/ ~ 4*(rA) J/) , 

(22) 



with the parameters w h(A) = = z rh{rA) = * h(M) = and 

v p E tE* , ; e* ; e 
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If li- 




Ce) (f) 

Figure 2: One loop diagrams contribute to r — > l\ decay due to the neutral Higgs bosons h 
and Aq in the 2HDM. % represents the internal lepton, /j~ (Z^~) outgoing lepton (anti lepton), 
internal dashed line the ho and A fields. 
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Figure 3: m r dependence of the BR (r — > r ± / u ± ). The solid-dashed lines represent the 
BR( r -> f±f) for £* >TT = 100 GeV, £* T „ = 10 GeV- £* >TT = 10 GeV, ^ = I GeV. 
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Figure 4: m r dependence of the BR 
) for = 100 GeV, & = 



r± e ± 



r -> If If 
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The solid-dashed lines represent the BR(r 



dashed line represents the BR 
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10 GeV, £f re = 0.1 GeV. The small 



= !GeV,a Te = 0.1GeV. 
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Figure 5: £§ TT dependence of the BR(r —r, ^ j ±u± <^ N:Tfi 
dashed lines represent the BR for the radion masses m r = 200 — 500 — 1000 GeV. 
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Figure 6: £^ TT dependence of the BR(r — > r^" for £ 



N,re 



0.1 GeV. The solid-dashed-small 



dashed lines represent the BR for the radion masses m r = 200 — 500 — 1000 GeV. 
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